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Abstract 



The distribution of transverse energy, Et, which accompanies deep-inelastic 
electron-proton scattering at small x, is predicted in the central region away 
from the current jet and proton remnants. We use BFKL dynamics, which 
arises from the summation of multiple gluon emissions at small x, to derive 
an analytic expression for the Et flow. One interesting feature is an x~'' in- 
crease of the Et distribution with decreasing x, where e = (3as/7r)21og 2. We 
perform a numerical study to examine the possibility of using characteristics 
of the Et distribution as a means of identifying BFKL dynamics at HERA. 

PACS numbers : 13.60.Hb, 12.38.Bx, 12.38.Qk 
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I. INTRODUCTION 

The experiments at HERA measure deep-inelastic electron-proton scattering in the previ- 
ously unexplored small x regime, x ^ 10~^. As usual, x is the Bjorken variable x = Q^/2p-q 
and = — g^, where p and q are the four momenta of the incoming proton and the virtual 
photon probe respectively. The deep-inelastic observables reflect the small x behaviour of 
the gluon, which is by far the dominant parton in this kinematic region. In particular the 
small X behaviour of the structure function F2{x, Q^) is driven by the g ^ qq transition. At 
such values of x soft gluon emission and the associated virtual gluon corrections give rise 
to powers of log(l/x), which have to be resummed. To leading aslog(l/x) order, this is 
carried out by the BFKL (or Lipatov) equation [jl| which, in a physical gauge, effectively 
amounts to the summation of gluon ladder diagrams. Two characteristic features of the 
result of this procedure are, first, the a;"'*' growth of the gluon distribution xg{x, Q'^) as x 
decreases, with A ~ 0.5, and, second, the relaxation of the strong-ordering of the transverse 
momenta, k^, of the gluons along the ladder (which is characteristic of the "large" x, large 
Q2 "Altarelli-Parisi" regime j . 

The recent measurements of F2{x,Q'^) at HERA show a rise with decreasing x 
which is entirely consistent with the x~^ BFKL perturbative QCD prediction Q. However 
the observed small x behaviour of F2 can also be mimicked by conventional dynamics based 
on the Altarelli-Parisi equation [^]. Because of the inclusive nature of F2, it is not a sensitive 
discriminator between BFKL and conventional dynamics. For this purpose it is necessary 
to look into properties of the final state. Here we investigate whether or not the emitted 
transverse energy. Ft, which accompanies a deep-inelastic event (see Fig. 1), can be used to 
identify the BFKL dynamics. Due to the relaxation of the strong-ordering of the gluon fc-p's 
along the chain we expect to find more transverse energy in the central region (between the 
current jet and the proton remnants) than would result from conventional evolution. Indeed 
a Monte Carlo study, in the context of heavy quark production, of the Ft distribution as 
a function of rapidity hints at an increase of the Ft flow if QCD "all-loop" dynamics 
(which incorporates BFKL effects) is employed, see Fig. 9 of Ref . . 

Now the Ft emitted by the gluon chain is given by the integral over the inclusive dis- 
tribution for the emission of one gluon weighted by its Ft- Thus, for a given element of 
(x, Q^), the distribution of transverse energy as a function of loga;^ is given by 

where a denotes the cross section in the (x, Q"^) interval, and where xj and kj are the 
longitudinal momentum fraction and the transverse momentum carried by the gluon jet in 
the photon-proton centre-of-mass frame. The variable logXj is closely related to rapidity, 
see Sec. |IIB| . In terms of the proton structure functions, the differential cross section in Eq. 
is 

da Ana"^ 
dxjdk'j xQ^ 

where y = Q'^/xs with y/s being the centre-of-mass energy of the colliding ep system. 



dF2 1 2^(2x^1 



dx^dk'i ~^ 2^ dx^dk'i 



:i.2) 



2 



With these definitions, we can calculate the energy fiow, XjOEt / dxj of Eq. ( |1.1| ), by using 
the diagrams of Fig. 2 to evaluate the differential structure functions. We have []To|-|T^ 

dxjdk] ^ / "ttA^ / ("tt If^) kl)6^'^\kj - - kp) (1.3) 

with i = 1,2, where Ti and / describe, respectively the soft gluon resummation above and 
below the emitted gluon of Fig. 2. That is jFj and / satisfy BFKL equations, provided of 
course that x/xj and Xj are sufficiently small. The BFKL equation for / may be written in 
the differential form 



\k'^-k^\ (4A;'4 + fc4)i_ 



K®f, (1.4) 



with a similar form |T3[ for jFj. As usual, / is the unintegrated gluon distribution in the 
proton, that is with the final integration over dk^ not performed. The expression in brackets 
in ( |1.3D arises from the (square of the) BFKL vertex for real gluon emission, where kp, k^ 
and kj are the gluon transverse momenta in the photon-proton frame. 

We may use Eqs (|l.l|)-(pr^), together with the solutions of the BFKL equations for / 
and J-'i, to determine the Et fiow, OEt / dlogXj, which accompanies deep-inelastic scattering, 
provided x -C a:j -C 1. Valuable insight can be obtained by using the analytic solutions of the 
BFKL equation for fixed a^. Therefore, before presenting the full numerical calculation, we 
perform in Sec. |I| an analytic study to establish general characteristics of the Et distribution. 
First, we make the simplifying assumption of strong-ordering at the vertex for real gluon 
emission; in fact in addition to the "normal" strong-ordering, fc| S> fcp, in the expression 
in brackets in Eq. ( |1.3| ), we also need to include the contribution from the "anomalous" 
ordering k^ ^ kj. Then we present an exact analytic study by using the BFKL form of the 
vertex, which incorporates all orderings automatically. In Sec. |TV| we present the predictions 
of the Et distribution obtained when the BFKL equation is solved numerically. At each 
stage we compare with the Et fiow pattern predicted by coventional dynamics based on the 
Altarelli-Parisi equation. Section |V] contains our conclusions. 



II. CHARACTERISTICS OF Et FLOW: AN ANALYTIC STUDY 

At first sight the prediction for the shape of the Et distribution as a function of log xj 
looks simple. The leading small z behaviour of the solutions of the BFKL equation for fixed 
ttg have the form 

f{z,k'^) z-^, with A = 4a, log 2 (2.1) 

where a, = 3as/7r, and similarly for jFj. Thus for x/xj and Xj sufficiently small, we see that 
the relevant combination which occurs in Eq. ( |1.3| ) has the behaviour 

-^,(x/x„ ...)/(x„ ..) ~ {-] (x,)-' (2.2) 
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which appears to give an Et distribution which is independent of xj, and also of x. Thus well 
away from both the current jet and the proton remnants, we anticipate a plateau distribution 
of Et as a function of \ogXj, which remains unchanged as x varies. 

However the situation is more subtle. We have overlooked the x and Xj dependences 
which arise from the integrations over the transverse momenta. To see this we need to look 
more closely at the analytic solution of the BFKL equation. Recall that the solution may 
be obtained by taking the Mellin transform / of /: 

f{z,e) = — / {erf{z,u;)du; (2.3) 

f{z,uj) = / ie)-^-'fiz,e)de (2.4) 
Jo 

with c = |. If we substitute Eq. ( p^) into ( |1.4|) we obtain 

= Kico)f{z,co) (2.5) 

where K, the transform of the BFKL kernel, is 

K{uj) = a,[2^(l)-^(cj)-^(l-cj)] (2.6) 
with \l'(u;) = r'(ci;)/r(u;). On substituting the solution of Eq. (^^ ) into Eq. (^^ ) we have 

fiz,^) = — I duj{erf{l,uj)z-^^-\ (2.7) 



2m 



c—too 



As uj varies along the contour of integration, the maximum value of K{uj) is found to occur 
at C(j = i. The region lu; ~ ^ therefore dominates the small z behaviour of /. To obtain the 
leading behaviour we expand the various terms in Eq. ( |2.7| ) about the point = | and find 



.2,^^-x(k'\^ foil) ..J-^og'ikypy 



f(z, e) ^ 7^ r exp "° ^ ; ^ / , (2.8) 

^ ' [27rA"log(l/^)]i 2A"log(l/^) y ' ^ ' 

with an analogous relation for the J^i 

where /o(|), ^i''°'*(|), ^^"^ ^ given in terms of the boundary conditions. The quan- 

tities A and A" arise from the expansion of KioS) of Eq. (|2l6| ) about the point u = \: 

k{uj) = \ + lM\u-\f + ... (2.10) 

with A given by Eq. ( |2.1| ) and A" = 0:^28^(3), where the Riemann zeta function ^(3) ~ 1.2. 

Note that the scale of the JFj is provided by and also that, with our definition ||13|| , 
the have the dimensions of Q^. In Eq. ( |2.8|) we have used a slightly different notation for 
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the function f{l,uj) at the boundary (which for simphcity we take to be 2; = 1). We have 
introduced a dimensionless boundary function /q defined by 

/(i,^) = (kir-Moo). (2.11) 

Similarly for the JFj we have introduced dimensionless boundary functions ^j;^^ defined by 

Mhu,Q') = iQ'Y~-ft\co). (2.12) 

To estimate the characteristics of the Et flow we may omit the longitudinal structure 
function and assume 2xFi = F2. It is straightforward to include the small correction arising 
from Fl = 2xFi — F2. Thus (|1 . 1|) becomes 



dEr 1 /• „2 dF 



where the integrand is specified by ( |1.3| ) with / and JF2 given by Eq. ( |2.8| ) and Eq. ( p.9| ) 
respectively, for sufficiently small Xj and x/xj. Thus to determine the Et flow we must 
perform the integrations dk"^, cPkp, cPk^ over the transverse momenta, subject of course to 
the conservation of transverse momentum as embodied in the delta function in (11.31). 



A. Simplified analytic treatment 



We can obtain valuable insight by making the assumption that the transverse momenta 
are strongly-ordered at the vertex for real gluon emission. At first sight, we would expect 
that this meant just keeping the contribution with /c| ^ k^ in Fig. 2. However for the E^ 
weighted distribution we shall find that there is an equally big contribution coming from the 
region with "anomalous" ordering k^ ^ k"^. We evaluate these two contributions in turn. 

For the "normal" strong-ordered case with fc? ^ k^ we have k"^ ^ /c|; then Eq. ( |1.3|) 
simplifies, since the cPkp/kp integration over / just gives the integrated gluon distribution 
Xjg{xj, kj). In fact in this approximation we recover the formula for deep-inelastic + ener- 
getic jet events 



(2.14) 



see, for example, Eq. ( p.4|) of Ref. [|T^]. Using the explicit solutions (|2.^ ) and (|2.9| ), and 
inserting ( ^^.141 ) into ( |2.13| ), we obtain 



dE, 



X 



(a) 
T 



^ dxj 



2a 



s -X 
-X, 



X 



-A 



Q'k 



2i.2\ 2 



_ /o(|)-^f (|)/(x,x,,g^) 

Co ) [27rA"log(l/xj)27rA"log(a;j/a;)]^ 



(2.15) 



where 



/ log^(A:|/P) log^(A:|/ag^) 



°° dk"^, Ik; 



kj (P)^ 



exp 



2A"log(l/xj) 2A"log(xj/x) 



(2.16) 
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The factor 2 in Eq. ( p.l5|) arises when we translate formula ( |2.8| ) for f{xj, k^) into one for 
Xjg{xj, kj). This factor can be inferred by using Mellin transform techniques and noting 
that the main contribution occurs when l/u ^ 2. 

The other contribution, which we denote XjdEj^^ /dxj, coming from the "anomalous" 
ordered region k^ ^ k^ can be readily shown to be equal to xjOe!^'' /dxj. To see this, we 
note that for this case k^ ~ kp and so we have 

= j^,2{kl)'^:F,{x/x,^lQ-)f{x,kll 



where the factor 2(kp)^ arises from 



j\k]rUk] = 2{ki] 
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If we insert Eqs ( |2.8|) and (|2.9| ) then we find XjdE^^ /dxj is also given by Eq. ( |2.15|) . This 



so-called anomalous contribution |T^ comes from ordering of transverse momenta that are 
opposite to that in the ordinary Altarelli-Parisi evolution of the structure function. We 
have found that in the small x limit and for fixed (albeit large) Q^, that it gives a con- 
tribution to the Et flow equal to that arising from the normal strong-ordering. Why does 
this "anomalous" contribution not occur, for instance, in the conventional double leading 
logarithm (DLL) limit of, say, F2{x, Q^)? In the total cross section (or in an inclusive struc- 
ture function) the dominant or leading "anomalous" contribution is cancelled by the virtual 
corrections which lead to gluon reggeisation p]JTl[| (or, equivalently, to the "non-Sudakov 
form factor" of Ref. |0-§])- The cancellation ensures that, at small x, we do not encounter 
the double logarithmic terms of the form aslog^(l/x). Moreover the virtual corrections 
guarantee that the conventional DLL terms log(l/x) logQ^ only come from the usual 
strongly-ordered region of transverse momenta, i.e. /Cp -C -C Q^. It may also be inter- 
esting to observe that, unlike the unweighted distributions, the Et flow is an infrared finite 
quantity since the factor of |kj| in the integrand of ( |1 . l\ j removes the unintegrable singu- 
larity which would otherwise occur at A;| = 0, see ( p..3|) . After this aside, we now return to 
complete the analytic calculation of the Et flow. 

The integral / of Eq. (|2.16|) can readily be evaluated if we change the variable from kj 
to t = \og{kj/k'^). Then Eq. (|2.16|) becomes 



— oo 



dt exp f- - - - -JLJl\ (2.17) 

^V2 2A"log(l/a;j) 2X" log{xj/x) J ^ ' 



where T = log(a(5^/fc^). If we evaluate this Gaussian integral, and insert the result into 
Eq. ( p^.l5D and into the identical expression for Xj Oe!^^ /dxj, we finally obtain the following 
analytic prediction for the Et flow 

OEt dEl^\ dE? .-,^2p,i-. f e[log{x]/x)+AT/X'f \ ,^ 

provided x/xj and sufficiently small, where e = A"/32. 
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To obtain Eq. ( p.l8|) we have used the BFKL prediction for F2 in the leading log 1/x 
approximation 0] 

= J^fix,e):Fi'\e,Q') (2.19) 

where jFg^"* is the quark box (and crossed box) contribution shown in Fig. 2(b), and where 
the integration over the gluon longitudinal momentum has been performed. If we insert 
expression (|2.8|) for / and we take the small x limit (where the Gaussian factor may be 
neglected), then the integral becomes proportional to the Mellin transform ^2^^ at 00 = |. 
Using Eq. ( |2.12| ) we obtain 

F, = ( %] ' , . (2.20) 

\kij [27rA"log(l/a;)]^ 



This result for F2 has been inserted into Eq. ( p.l5| ). 

Equation ( |2.18| ) for the Et flow applies in the central region away from the current 



jet and the proton remnants, that is when Xj/x ^ 1 and l/xj ^ 1. Contrary to our 
initial expectation, (|2.2| ), the expression (|2.18|) depends on x, and also on xj. The main 
characteristics of the flow can be readily identified from Eq. ( p.l8| ). We see that the E^ 
distribution has a broad Gaussian shape in logXj (or in rapidity, see Sec. II B| below) with 



(i) a peak which moves away from xj ~ 0:2 as increases. 



(ii) a width which grows as J\og{l/x) decreases, 

(iii) an x~^ dependence such that it grows with decreasing x, 

(iv) a {Q'^)^ dependence such that it grows with increasing Q^. 

Although the shape of the distribution is specified, the normalisation is controlled by log(A;^), 
the parameter which gives the centre of the input (non-perturbative) log /e^ Gaussian distri- 
bution of /, see Eq. (O) and Ref. §. 



B. Et flow as a function of rapidity 

We have calculated the Et flow as a function of log xj, a variable which is closely related 
to the rapidity, y. The result, OEt / dlogXj, is the Gaussian form of Eq. ( |2.18|) . For practical 
purposes the rapidity distribution dET/dy is more useful, so here we show how to evaluate 
the integral I of Eq. ( p.l7[ ) in terms of y, rather than logXj. In this way we will find that 
the Et flow as a function of rapidity also has a Gaussian form. 

We define the rapidity y in the photon-proton cm. frame in the standard way 

y^|log(i|). (2,21) 

which leads to the following relation between y and log Xj 
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logx, - i/-i(T-t)-ilog(l/x). (2.22) 

In terms of y, the integral / of Eq. ( p.l7| ) assumes the (non- Gaussian) form 

'•fe ft _ e _ it-Tf \ ( \ 

u ^ """"P 1,2 " \"[-2y + T-t + log(l/x)] ' \"[2y-T + t + log(l/x)] ) ^^"^^^ 

where the hmits 

ti,2 = T-2yTlog(l/a;). (2.24) 
Now, to obtain an analytic form we approximate the t factors in the denominators in Eq. 

i = 4e log(l/x) + iT, 

which denotes the value of t for which the integrand in Eq. ( |2.17| ) is a maximum in the 
leading log(l/x) approximation. It is straightforward to evaluate the integral approximated 
in this way. We find 

dy V log(l/a;) / 

Thus we obtain a Gaussian Et distribution as a function of rapidity y, identical in structure 
to that in terms of logXj, except for the two non- leading terms in the exponential. 



C. Improved analytic treatment 



Equation ( p.l8| ), though valuable to see the origin of the x and xj dependence of Et, 
is based on assumptions which are not necessary even for an analytic treatment. In this 
subsection we show that it is not necessary to assume (i) strong-ordering at the gluon 
emission vertex and (ii) the approximate Gaussian diffusion patterns of Eq. ( p.8| ) and Eq. 
( p.9|) . The procedure is to evaluate (|1.1| )- (|1.3| ) using Mellin transform techniques. We have 



OEt 
dxj 



1 



1.2 



a. 



COS( 



1 



dkl 



-H{k'^/k'p)J^,{x/x,, fc^, Q')f{x„ k'^). 



^1 J ^ ik'.k',) 

The function H, introduced as above, is now replaced by the inverse Mellin transform 



(2.26) 



2m 



duH{u)U 

-IOC \ ftp ^ 



(2.27) 



with c in the range — | < c < Proceeding in this way we see that the kp and k^ 
integrations, over / and JF2 respectively, now factorize and, moreover, are simply of the form 
of the Mellin transforms / and defined as in Eq. (|2.4| ). Thus Eq. ( p.26|) becomes 
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bet 1 1 

OXj r2 



i i fC+lOO _ _ 

- F 2^7™ d''H{u)J^,{x/x,,\-u,Q')f{x,,l + u). (2.28) 



The Mellin transform f{z,u) satisfies the transformed BFKL equation ( p.5| ) which has the 
solution 

fiz,u) = (A;2)-/oHz-^(-), (2.29) 

see Eq. ( p.ll[ ). Similarly the solution of the equation for JF2 is 

:F2{z,u;,Q') = m'--fi'\u)z~^^^\ (2.30) 

recall Eq. ( ^^.121 ). The explicit form of the function ^2'\uj) can be simply deduced from, 
for example, Appendix B of Ref. Inserting solutions Eq. ( |2.29|) and Eq. ( |2.3(]|) into Eq. 



( [2.2^ ), we obtain 

= IT- duH{u) ,,,,, „ J^2 (l-'^)/o(i + '^) — 



OXj F2 2mJc-ioo (^0)4+^^ 



(2.31) 



We see that in the limit where Xj/x 3> 1 and 1/xj 3> 1, and in particular in the central 
region away from the current jet and the proton remnants where Xj ~ x^ (or to be more 
precise where log(l/a;j) ~ log(l/x2)), that the integral is dominated by contributions from 
the region with i/ ~ 0. We therefore expand the relevant quantities about i/ = and evaluate 
the resulting Gaussian integral. We find the master analytic formula 

rfn^P^k ( {Wosjxy x)+fr \ 

Xj— — = CiQ k^j^x exp ^ , (2.32) 

'dx, ' 2A"log(l/x) r 



which may be compared with result ( p.l8| ) obtained by the simplified analytic treatment. 
Here T = log{Q'^ /k'^). We see we have recovered the same general structure as before, that 
is a Gaussian Et distribution as a function of logXj, with the characteristics (i) to (iv) listed 



in Sec. P A| . The quantities are now more precisely determined. First the parameter e, which 



controls the growth of the distribution with decreasing x, is given by 

6 = Ka)-ka) = 2«,log2, (2.33) 



which would reduce to the previous result, e = A"/32, if the quadratic form Eq. ( p.lO|) of 



K{uj) were valid to u; = |. The parameters A' and A" denote the derivatives of iT at a; 

A' = K'i\), X" ^ K"i\), (2.34) 
and finally the normalization 
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where 



■^4/1^ 



^(0) 



27rP(i) 



8.75a,. 



(2.36) 



As before we can translate the Gaussian Et distribution in logxj to one in terms of 
rapidity y. Again the distribution turns out to have the form of Eq. ( p. 321) , 



-{\X'\2y + f + fX'Hog{l/x) - l\\'\fy 



dy - V- / ^ y 2A"log(l/x) 
but we see that there are additional, non-leading, terms in the exponential. 



(2.37) 



III. Et flow with STRONG Kt ORDERING IMPOSED 



We compare the BFKL forms of the Et distribution, Eq. ( p.l8|) and Eq. ( p.32|) , with 
that obtained from conventional dynamics based on Altarelli-Parisi evolution in the small 
X limit. To make a strict comparison we take fixed a<j. Thus rather than inserting into Eq. 
( p.l4| ) the BFKL solutions (^78|) and (|2.9|) , we replace them by the double leading logarithm 
(DLL) forms 



exp 



(3.1) 



k"^ exp 2 



(3.2) 



We substitute these forms into Eq. ( p.l4| ) and Eq. ( |2.13| ), and change the fc| integration 
variable to t = log(/c|/(5o)- then obtain 



QEt as , 2^ 
where T = log(QVQg 



2 f cit exp I - + 2 
\ 2 



as{T -t) log 



X 



+ 2 



Ust log ■ 



Xj 



It is informative to evaluate the Et distribution for that is for 

log A = log- = 



Xa 



X 



(3.3) 



(3.4) 



a point "midway" between the current jet and the proton remnants. This choice facilitates 
an analytic study and allows us to gain insight into the x and dependence of the Et flow. 
When \og{xj/x) = log(l/xj) we can evaluate Eq. ( ^73|) using saddle point methods. We find 



dEj 



dxi 



as{Qlf'expnf3T+^2asTY + - 



(3.5) 
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where (3 is a known function of the ratio T/2asY. To be precise 



P 



1-4 



l + Jl + T/2a,Y 



-21 2 



(3.6) 



The — in the square brackets in Eq. ( |3.5| ) arises from substituting 



exp 2 



(3.7) 



into Eq. ( |3.3|) , rather than using Eq. (|2.20| ). After some algebra, it is possible to show that 
the approximate estimate of the integral in Eq. ( |3.3| ) predicts a slow increase of Et with 
decreasing Xj. 

We inspect the properties of the Et flow, Eq. ( p.5| ), in two different limits. First in the 
double leading logarithm limit, where agTY ^ 1 but a^T ~ 1 and agY ~ 1, we see that 
/3 — i> 1, and so 



^E^ 



exp 



2-l)J2a,TF . 



(3.8) 



Thus in the DLL limit we find that the Et distribution grows essentially as ((5^)2 with 
increasing . In the second limit, l/x — *• oo at finite Q^, we see that ^5 ^ \. Thus we find 



dEn 



(3.9) 



Moreover it is possible to show from Eq. ( |3.5| ) that with decreasing x the distribution de- 
creases towards its final limit, which is indicated by Eq. ( p.9|) (modulo slowly varying loga- 
rithmic factors). This should be contrasted with the BFKL behaviour of Eq. ( |2.32| ) which 
increases as x"*^ with decreasing x. 



IV. CALCULATION OF Et FLOW 

Above we have obtained analytic forms for the leading behaviour of the Et distribution 
in the Xj 0, x/xj limit, for fixed a^. These expose characteristic features of the 
expected Et flow in deep-inelastic scattering in the central region between, but well away 
from, the current jet and proton remnants. To obtain a more realistic estimate we solve the 
appropriate BFKL equations numerically, and use running a^. From Eqs ( |1.3| ) and (|2.13| ) 
we see that the Et flow, for a given x, Q^, can be expressed as 




(4.1) 

In this paper we focus on the central region with Xj < 10~^ and x/xj < 10~^. Thus to 
calculate the Et flow we have to solve three different BFKL equations, one for / which 
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resums the lower gluon ladders of Fig. 2(b), one for T2 which resums the upper gluon 
ladders, and finally one for F2(x,(5^). To be precise we determine f{xj,kp) for Xj < 10~^ 



by solving the BFKL equation, as described in Ref. |T3[, starting from a gluon distribution 
at Xj = 10~^ obtained from the MRS set of partons [0. The quantity ^^2/^ may be 
identified with the structure function of a gluon of (approximate) virtuality k"^, integrated 
over its longitudinal momentum. The function ^-2(2;, /c^, Q^) is determined for z < 10~^ by 
solving the BFKL equation, as described in Ref. , starting from the boundary condition 
^2 = ^2°^ at 2; = 10~^, where ^-"2^^ is the contribution resulting just from the quark box 
(and crossed box) in Fig. 2(b). Finally F2{x,Q'^) is determined by the convolution of the 
gluon ladder with the quark box (and crossed box), which symbolically may be written 
F2 = / (g) as described in §. 

The infrared cut-off on the transverse momenta integrations is taken to he = 1 GeV^ 
throughout, a value which makes the calculated values of F2 consistent with the recent 
observations at HERA 0,^. The predictions for the Et flow are less sensitive to the choice 
of the cut-off than those of F2 . The ultraviolet cut-off on the integrations over the transverse 
momenta is taken to be 10"^ GeV^; the results for /, J^2 and F2 are not sensitive to reasonable 
variations about this value. We shall see below, however, that the E'j^-weighted integration 
of Eq. ( ^4.1| ) is sensitive to the ultraviolet cut-off as xj decreases towards the current jet. 

Numerical results for the Et flow are presented in Fig. 3 for = 10 GeV^ and values 
of X = 10"*^, 10"'^, ...10~^. We show results for such low values of x so as to establish a 
sufficiently large central region in order to compare the general features of the distribution 
with the expectations of the leading order analytic treatment of Sec. |I^ From Fig. 3 we 
indeed see a broad Gaussian shape in logXj, with a peak below and a width which 

grows with decreasing x. Moreover the full numerical treatment shows an increase of the 
height of the distribution with decreasing x, but not as rapid as anticipated by the leading 
order analytic result. 

We also show in Fig. 3 the Et distributions obtained using Altarelli-Parisi (AP) evolution 
for two values of x. These confirm the expectation that they decrease with decreasing x. 
We come back to the BFKL comparison with AP when discussing Fig. 5, but first in Fig. 4 
we explore the sensitivity to the choice of the ultraviolet cut-off in Eq. (|4.1| ). 

Having established the general features of the Et flow, we focus on the distribution for 
X = 10~^, a value which is accessible at HERA. In Fig. 4 the sequence of dashed curves show 
the effect of choosing the ultraviolet cut-off in Eq. { ^A\j to be /c^^ = 10^, 10^, ...10"^ GeV^, and 
the continuous curve is the result obtained if the cut-off is chosen to be Q'^/z as implied by 



energy-momentum conservation |17], where z = x/xj. At relatively low values of Q^, such 
as = 10 GeV^, which are relevant to experiments at HERA, energy conservation has a 
significant effect on the Et flow arising from BFKL gluon emission. Fig. 5 compares the Et 
flow for g2 = 10 GeV^ with that for = 100 GeV^ with, in each case, kl^ = Q^/z. Also 
shown are the expectations from Altarelli-Parisi evolution. As anticipated we see a significant 
difference between the Et flow arising from BFKL and AP gluon emission, particularly at 
lower values. 

We may use Eq. ( |2.21| ) to translate the Et flow as a function Xj into a distribution in 
terms of the rapidity y in the virtual-photon-proton cm. frame. For a reasonably flat Et 
distribution a good estimate of y is obtained if we insert the average value of k^ into Eq. 
( p.21| ). For = 10 GeV^ we see that Et ~ 2 GeV per unit of rapidity and so we take 
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kj=4: GeV^. The result for y is the scale shown above Fig. 5. If the Et distribution 
observed by the HI collaboration [|l^ is compared with Monte Carlo expectations, then 
reasonable agreement is found for x > 10~^, but for x < 10~^ more Et appears to be 
present in the region between the current jet and the proton remnants than is implied by 
the favoured Monte Carlo calculations (which at present do not allow for BFKL emissions). 
The deficiency is most evident in the rapidity region xj ^ 10~^, which extends beyond the 
central region shown in Fig. 5. However it is encouraging that a simple extrapolation of the 
BFKL curves of Fig. 5 into this region indicates a sizeable enhancement of the predicted E^ 
distribution over that expected from conventional AP dynamics. 



V. CONCLUSIONS 

Our main objective is to find characteristics which are unique to BFKL dynamics at small 
X. One striking feature of BFKL dynamics is the very strong a;"'^, with A ~ 0.5, increase 
of the deep-inelastic structure function E2{x, Q"^) as x decreases in the HERA regime 0,|19 



However this increase can be mimicked by conventional Altarelli-Parisi evolution from a 
very low starting scale Ql = 0.3 GeV^ [|^. Global or inclusive quantities like structure 
functions do not provide conclusive tests of BFKL dynamics and it has been advocated 
that the BFKL effect can, in principle, be more definitively identified by studies of the final 
state in deep-inelastic scattering. One possibility is Mueller's proposal [^] of measuring 



the x/xj dependence of deep- inelastic events containing a forward-going identified jet of 
longitudinal momentum Xjp. Another, which we have studied here, is to measure the Et 
flow accompanying deep-inelastic events at small x. A distinctive feature of BFKL dynamics 
is the diffusion pattern of the transverse momentum flow, in contrast to the strong-ordering 
of kT towards which is characteristic of conventional Altarelli-Parisi dynamics. Indeed it 
is the absence of strong /cT-ordering which leads to a strongly enhanced Et distribution at 
small X. We analysed this effect both analytically and numerically, focussing attention on 
the central region between the current jet and the proton remnants. 

In Sec. |Tl|we performed a detailed analytic analysis of the Et flow. We found that BFKL 
dynamics with fixed Og gives a broad Gaussian Et distribution as a function of rapidity (or 
rather logxj) which grows as x~'^ with decreasing x where e = (3as/7r)21og2. This should be 
contrasted with the much smaller Et flow obtained assuming strong fcr-ordering, which gives 
an Et distribution that decreases with decreasing x, for fixed Q^. In Section 4 we confirmed 
these qualitative features using numerical solutions of the appropriate BFKL equations. As 
expected we found that the Et flow which accompanies deep-inelastic events is much larger 
for BFKL than for conventional AP dynamics. We also showed, at the low values relevant 
to HERA, that energy conservation should be used to limit the energy-weighted integrations 
over transverse momentum. 

Taking, for example, deep- inelastic events with x = 10~^ and = 10 GeV^ we found 
that the Et flow in the central rapidity region {xj < 10~^, x/xj < 10~^) between the 
current jet and the proton remnants is about 2 GeV per unit of rapidity. This result is 
very encouraging for the identification of BFKL dynamics. However for a quantitative 
comparison with experiment we will need to include contributions to the Et flow coming 
from the hadronization of the produced partons, not only from the BFKL ladder, but also 
from the partons radiated from the current jet. Also there is a contribution radiated from 
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the colour string which connects the current jet to the proton remnants. An analysis of this 
problem is under study and its results will be reported elsewhere. 
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FIGURES 



FIG. 1. Deep- inelastic scattering at small x, accompanied by soft gluon emission with total 
transverse energy TjExi- 

FIG. 2. (a) Diagrammatic representation of Eq. (b) explicit display of the gluon ladders 

which are resummed by the BFKL equations for / and J-i. 



FIG. 3. The Et flow as a function of logxj, calculated using BFKL forms in Eq. ([4.1|), for 
X = 10~®, 10"'', ...10~^ and = 10 GeV^, in the central region {xj < 10~^, x/xj < 10~^) between 
the current jet {xj ~ x) and the proton remnants {xj ~ 1). For comparison, the dashed curves 
show the Et flow calculated using Altarelli-Parisi evolution for x = 10^^ and 10~^. 



FIG. 4. The sensitivity of the Et flow (for x = and = 10 GeV^) to the choice of the 
ultraviolet cut-off k^^ imposed on the integrations over the transverse momenta in Eq. ( |4.1| ). The 
dashed curves correspond to k'^-^ = W'^ , 10^, ...10^ GeV^, and the continuous curve to k^^ = Q'^ / z 
with z = x/xj. 



FIG. 5. The Et flow for x = 10"^ calculated using BFKL forms (continuous curves) and 
conventional AP dynamics (dashed curves) for two different values of Q^, namely = 10 and 
100 GeV^. The BFKL predictions are obtained from Eq. (4.1) with k'^^ = Q'^/z with z = x/xj. 
The variable logXj is related to the virtual-photon-proton cm. rapidity, y, via Eq. (2.21). The 
approximate rapidity scale for = 10 GeV'^ is shown above the figure. 
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